Towards the grain boundary phonon scattering problem: an evidence for a 

low-temperature crossover 
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The problem of phonon scattering by grain boundaries is studied within the wedge disclination 
dipole (WDD) model. It is shown that a specific (/-dependence of the phonon mean free path for 
biaxial WDD results in a low-temperature crossover of the thermal conductivity, k. The obtained 
results allow to explain the experimentally observed deviation of k from a T s dependence below 
Q.IK in LiF and NaCl. 
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The effect of low-angle grain boundaries on the ther- 
mal conductivity, k, of LiF and NaCl over the temper- 
ature range 0.08-5 K has been investigated in Jl],^]. The 
main conclusions are the following: (i) the boundaries 
are sessile, (ii) the dominant phonon-scattering process 
comes from static strain fields caused by boundaries, and 
(iii) the experimental results are compatible with predic- 
tions of the theoretical model J3j where a grain bound- 
ary is represented as a wall of edge dislocations. In ex- 
periments, however, in addition to the expected behav- 
ior kT~ 3 — const, a remarkable increase in kT~ 3 below 
T* ~ 0.1K was detected. A similar deviation, but be- 
ginning near 2K, was observed in sapphire There 
is still no satisfactory explanation of this phenomenon. 
In particular, in jjj it was supposed that the measured 
increase can be caused by the onset of partial specular 
reflection from the lightly sandblasted walls. Similarly, it 
was suggested that [Q " a frequency independent scatter- 
ing mechanism should be present in these samples which 
becomes ineffective below IK" . 

In 1955 Klemens j| studied the problem of the scat- 
tering of lattice waves by grain boundaries within the 
Born approximation. Considering the grain boundary as 
an array of edge dislocations lying in the plane of the 
boundary, he found that the phonon mean free path is 
frequency independent. Hence a T 3 dependence of the 
thermal conductivity at low temperatures was associated 
with the boundary scattering. While this finding explains 
well the experimental results above some charac- 

teristic temperature, T*, it fails to describe the observed 
anomaly below T*. It is important to note in this con- 
nection that the result || was obtained under assump- 
tion that the dislocation wall is infinitely long. For a 
finite wall of well separated dislocations the problem of 
the phonon scattering becomes difficult and is still unre- 
solved. 

An alternative model for description of grain bound- 
aries has been presented in J^] . It was proposed that grain 



boundaries being rather rotational than translational de- 
fects can be described more naturally by disclinations. 
Moreover, what is important, the far strain fields caused 
by wedge disclination dipoles (WDD) were found to agree 
with those from finite walls of edge dislocations ||f7| . For 
this reason, the WDD-based model allows us to study im- 
portant effects due to finiteness of grain boundaries. 

Notice also that additional interest to this problem 
was inspired by recent consideration of disclinations 
and dipoles of disclinations in the context of metal 
glasses [||,[)| , graphite films [JlOj , and nanostructures [^l) . 
For example, an attractive model for a metallic glass pro- 
posed in Q results from disordering of a Frank-Kasper 
disclination network. Graphite films are expected to con- 
tain a number of disclination pairs due to the presence 
of the five- and the seven-membered rings ||lo|| . Discli- 
nations and disclination dipoles are of importance in 
graphite nanotubes which have attracted great interest 
recently. 

At the same time, it is known || that information 
about the nature of the principal imperfections in crystals 
can be extracted from a study of the thermal conductiv- 
ity. In particular, the temperature dependence of the 
thermal conductivity arising from defects is governed by 
the frequency dependence of their scattering cross section 
for lattice waves. 

In this Letter, we study the problem of phonon scatter- 
ing by grain boundaries in the framework of the WDD 
picture. These scatterers are shown to possess rather 
specific properties. In particular, the q-dependence of 
the phonon mean free path varies significantly for dif- 
ferent types of dipoles. We show that the experimental 
results |^J^] are in good agreement with the calculated 
thermal conductivity due to phonon scattering by biaxial 
WDD both above and below T* . 

Let us calculate a mean free path of phonons of fre- 
quency u> scattered by the potential associated with a 
static deformation of a lattice caused by straight WDD. 
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Following the generally accepted approach |||l2j we con- 
sider an effective perturbation energy due to the strain 
field caused by a single WDD in the form 



U(r) = hujjSpE AB , 



(1) 



where Tiuj is the phonon energy with the wavevector q, 
lu = qv s , v s is the sound velocity (it is assumed that three 
acoustic branches are equivalent), 7 is the Griineisen con- 
stant, and Eab is the strain tensor due to WDD. 

To simplify the problem, we assume that incident 
phonons are normal to disclination lines. The suitable 
geometry is chosen: disclination lines are directed along 
the z-axis, the dipole's arm is oriented along the x-axis. 
In this case, the strain matrix for WDD is known (see, 
e -g-j @)> an d U(r) takes the form: 



U(x,y)=B 



+ fa 



L 



1 (x + L) 2 +y 2 

2 11 (x - L) 2 + y 2 1 (x + L) 2 + y 2 



L 



(x - L) 2 + y 2 



(2) 



where B = hqv s ^fv(l — 2er)/(l — a), v is the Frank index, 
a is the Poisson constant, 2L is the dipole separation, 
and parameters l\ and fa specify displacements of axes of 
rotation from points where disclination lines pierce the xy 
plane. Notice that all possible types of WDD are included 
in Eq.(2). In general, it describes a biaxial WDD with 
arbitrarily arranged axes of rotation. For fa — fa = 2Z, 
fa = —fa, and fa = fa = one obtains the uniaxial WDD, 
the symmetrical uniaxial WDD, and the biaxial WDD 
with nonskew axes of rotation, respectively. 

According to Eq.(2), at chosen geometry the problem 
reduces to the two-dimensional case with the phonon 
mean free path given by 



A 9 1 = n, 



2tt 



(1 - cos e)^{e)de. 



(3) 



Here is an effective differential scattering radius, 

and rii is the areal density of defects. Within the Born 
approximation 5R(6*) is determined as [p"2[ 



qS 2 



2-kUv 2 



<q\U(r)\q> > 



(4) 



where all vectors are two-dimensional ones, S is a pro- 
jected area, the bar denotes an averaging procedure over 
a which defines an angle between p = q — q 1 and the 
x-axis. In other words, it means the averaging over ran- 
domly oriented dipoles in the xy plane. Evidently, the 
problem reduces to the estimation of the matrix element 
in Eq.(4) with the potential from Eq.(2). For this pur- 
pose, it is convenient to use the polar coordinates (r, (f>) 

U(p,a) = <q\U(f)W > 

= — J d 2 rexp{iprcos((f)~a)]U(r,(f)). (5) 



Here p = \p\ = 2q sin (0/2). Omitting the tedious calcu- 
lations (the details will be presented elsewhere) we write 
out the final result 



47rz 

U{p,a) = B — sin(pL cos a) 

p z 

H cos a coslpLcosa) 

P 



(6) 



where A; = h — fa- After averaging of \U(p,a)\ 2 over a 
and following integration in Eq.(3) with respect to 9 one 
obtains 



A" 1 = D 2 { V Lfn iq { z 2 []- + J 2 (2qL) 



z(z + ; 



(j 2 (2qL) + J 2 (2qL)) 



qL 



J Q (2qL)J 1 (2qL) \, 



(7) 



where D = irj(l — 2er)/(l — cr), z — Ai/L, and J n (t) are 
the Bessel functions. It should be emphasized that Eq.(7) 
is the exact result which allows to describe all types of 
WDD. Notice, that the behavior of A q in Eq.(7) is ac- 
tually governed by the only parameter 2L which charac- 
terizes the dipole separation. Before proceeding to im- 
portant applications of Eq.(7) let us consider two limits 
to early solved problems. As mentioned above, a biaxial 
WDD (more precisely for li = fa = 0) can be simulated 
by a finite wall of edge dislocations with Burgers vectors 
situated in parallel. As a consequence, with decreasing a 
dipole separation such dipole should be equivalent in its 
properties to a single edge dislocation with the Burgers 
vector b = AirLve y . Indeed, for small L one can obtain 
from Eq.(7) that A ~ g -1 . This is the well-known result 
for the phonon scattering by edge dislocations. In the 
opposite limit of large L, A q for a biaxial dipole becomes 
constant (see details below) in agreement with the result 
obtained for the infinite dislocation wall ||. 

In real materials, dipole separations which are of or- 
der of the grain size can lie in the broad range between 
50A (for nanocrystals) and 10 5 A (for polycrystals). As 
it follows from Eq.(7), two specific regimes of scatter- 
ing appear depending on the wavelength A of an incident 
phonon in comparison with 2L. As would be expected, 
the changing in behavior of A g occurs at A ~ 2L. Thus, 
for thermal phonons, we can estimate the typical tem- 
perature of transition 



T* 



2Lk f 



(8) 



where ks is the Boltzman constant. 

An interesting consequence of Eq.(7) is the conclusion 
that despite the seemingly similar nature of scatterers 
the behavior of A q differs remarkably for different types 
of dipoles. Let us examine two limiting cases. 
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1. Uniaxial dipoles, A; 
form 



2L (z = 2). Eq.(7) takes the 



A" 1 = 2D 2 (yL) 2 n l q(\ + J 2 (2gL) - jf (2gL) 



— J (2qL)J 1 (2qL) 
qL 



(9) 



According to Eq.(9), for gL > 1 one obtains A q ~ g -1 , 
whereas for <C 1 A g ~ g -5 . Thus, in the case of small 
wavelengths (in comparison with a dipole separation) the 
phonon scattering due to uniaxial WDD behaves like that 
for edge dislocations. In the limit A 3> L the scattering 
of phonons by uniaxial WDD is found to have a strong 
q -dependence, even stronger than the known Rayleigh 
scattering of phonons by point impurities. It should be 
noted that such behavior is compatible with an overall 
view of uniaxial WDD as a strongly screened system . 
It is worthwhile to mention that the low-angle uniaxial 
WDD can be simulated by a finite wall of edge disloca- 
tions complemented by two additional edge dislocations 
at both ends of the wall. The sign of these dislocations is 
opposite to that of dislocations in the wall and absolute 
values of Burgers vectors are equal to b = 2Ltan(7rz/). 
Notice that just these two dislocations provide a screen- 
ing. A significantly different picture arises in the case of 
the biaxial WDD with l\ = l 2 - 

2. Biaxial dipoles with A; = [z = 0). Eq.(7) trans- 
forms to 



A" 1 = 8D 2 (isL) 2 n iq (j 2 (2qL) + J x 2 (2gi) 



— Jo(2qL)J 1 (2qL) 



(10) 



In this case, in the long wavelength limit one gets A q ~ 
q^ 1 while for A < L we obtain that A q — ► const. It should 
be stressed that the appearance of the (/-independent re- 
gion distinguishes remarkably this scatterer from other 
types of WDD as well as from dislocations. 

Fig.l shows A q for three types of WDD at 2qnL = 
6 x 10 3 (or L w 2 x 10 3 A). One can clearly see the char- 
acteristic points where crossover occurs. In accordance 
with our estimation in Eq.(8), crossover takes place at 
q/q D ~ (1/6) x 1(T 3 or, respectively, at T* ~ 6/6000 
where is the Debye temperature. In the general case 
of arbitrarily arranged axes of rotation, < A < 2L, the 
mean free path interpolates between two limiting curves. 
The constant in Eq.(7) is chosen to be D 2 {vL) 2 niq£, — 




FIG. 1. Phonon mean free path Aq as a function of scaled 
frequency q/qD at 2quL — 6 x 10 3 for A; = (solid line), 
Ai = 2L (dashed line), and A; = 0.5L (dotted line). The 
1.35xl0~ 5 cm, v = 0.023, D = 2.6, 

,5, 



parameter set used is: L 



ni = 1.8 x 10 cm , and v s = 4.8 x 10 cm /sec. 

Let us discuss briefly the applicability of the Born ap- 
proximation to the considered problem. Notice that, as 
in the case of dislocations, the far field of the potential in 
Eq.(2) behaves like 1/r. Thus, all arguments in support 
of the validity of the Born approximation for the problem 
of phonon scattering by dislocations (see, e.g., [p^| , p^[ ) 
are appropriate here as well. Nevertheless, it is useful 
to perform a rough estimation. In our case, by anal- 
ogy with unscreened Coulomb potential [EJ, we obtain 
2qL <C (1 — 2cr)/[(l — a)uj]. Thus, one can conclude that 
the Born approximation should be valid for description 
of low-angle WDD and/or small dipole separation 2L. 
Otherwise, one has to restrict the region of admissible q. 
One can expect, however, that the qualitative behavior of 
A q will remain unchanged even beyond the applicability 
of the Born approximation. Notice also that when axes 
of WDD are orientated randomly, one has to perform an 
additional averaging. However, the known result for dis- 
locations |^,|l^| indicates that such averaging should lead 
only to a modification of the numerical factor in Eq.(7). 

Let us estimate a contribution to the thermal conduc- 
tivity caused by the phonon scattering due to WDD. For 
this purpose, one can use the known kinetic formula writ- 
ten in the dimcnsionlcss form 



k%T* 



2ir 2 fr 



3 v 2 



e/r 



l e x (e x - iy 2 A(x)dx, (11) 



where x — fiu/hsT, and the specific heat capacity is 
chosen in the standard Debye form. We have restricted 
ourselves to the thermal phonons with T < O. Then, for 
uniaxial WDD one obtains that k ~ T~ 2 at low temper- 
atures while k ~ T _1 for T — > O. 

Let us examine Eq.(ll) in detail for biaxial WDD with 
A q from Eq.(10). In accordance with the above analysis, 
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the thermal conductivity should exhibit a crossover from 
k ~ T 2 to k ~ T 3 with T increasing. Although the 
crossover temperature To differs from T* , it is still deter- 
mined by the value of 2L. 
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exactly calculated. We have shown that the thermal con- 
ductivity exhibits a crossover from k ~ T 2 to k ~ T 3 with 
T increasing. Thus, the biaxial WDD is a good candi- 
date for the specific scatterer proposed in p| . The results 
obtained allow us to explain the deviation of the thermal 
conductivity from a Independence below O.IK observed 
in LiF and NaCl. We expect that this crossover has a 
universal character and should be observed in various 
materials with the pronounced granular structure. 
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FIG. 2. Reduced thermal conductivity due to WDD scat- 
tering, k x T~ 3 vs temperature T, calculated according to 
Eq.(ll) with the same parameter set as in Fig.l. Measured 
points for the boundary-limited thermal conductivity in LiF 
(from Ref.[l]) are indicated by triangles. 

Fig. 2 shows the reduced thermal conductivity calcu- 
lated theoretically according to Eq.(ll). We would like 
to stress that, within our scenario, an increase in kT~ 3 
below T* arises from the phonon scattering by biaxial 
WDD with li = I2 and is quite universal. It should be 
observed in polycrystals with appropriate sizes of grains. 

We considered the experimental results for LiF E10|, 
NaCl 0, and sapphire [Q. Experimental data due to 
boundary-induced phonon scattering in LiF from O are 
shown in Fig. 2. As is seen, the data are in a good agree- 
ment with our results for appropriate choice of model 
parameters. Notice that for a qualitative comparison we 
need the only parameter which is a size of the grain. For 
NaCl the thermal conductivity was found to have a sim- 
ilar behavior but with larger magnitude [El. In sapphire, 
the measured crossover temperature is essentially higher. 
This can be explained by smaller grain sizes. It is inter- 
esting to note in this connection that small L can take 
place in glasses where the crossover temperature T ~ IK 
and even higher. Thus, the dependence k ~ T 2 can be 
extended up to these temperatures. A detailed study of 
physics of dielectric glasses on the basis of biaxial WDD 
will be presented in a separate paper. 

To conclude, the grain-boundary-induced phonon scat- 
tering has been studied within the WDD-based model. 
The proposed model is shown to take into account the 
finiteness of the boundary which is important for ther- 
mal phonon scattering. The phonon mean free path due 
to scattering by static wedge disclination dipoles has been 
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